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CN , Abstract 



We consider the Schrodinger operator with a periodic potential on quasi-lD models of 
armchair single- wall nanotubes. The spectrum of this operator consists of an absolutely 
continuous part (intervals separated by gaps) plus an infinite number of eigenvalues 
with infinite multiplicity. We describe all eigenfunctions with the same eigenvalue. 
We define a Lyapunov function, which is analytic on some Riemann surface. On each 
Q_|. sheet, the Lyapunov function has the same properties as in the scalar case, but it has 

branch points, which we call resonances. In example we show the existence of real and 
complex resonances for some specific potentials. 



(N 

s 



1 Introduction and main results 

7— I ■ 
> 

ON \ Consider the Schrodinger operator J$? = —A + "V q with a periodic potential "V q on so called 

armchair graph T N , N ^ 1. In order to describe the graph T N , we introduce the fundamental 
(T) ■ cell T = UjeNeTj C M. 2 , where Tj = {x = rj + tej,t G [0, 1]} are edges of length 1, 

l> 
O 

o 



N m = {1,2, ..,m}, ex = e 6 = Vs), e 2 = e 4 = (l,0), e 3 = -e 5 = ^(1, -V3) 



ri = (0, 0), ? 2 =? 5 =? 1 + e 1 , ? 3 =? 6 =7 2 + e 2 , r 4 = r 3 + e 3 . (1.1) 

We define the strip graph T N by 

f N = U (n , fc)eZxNiV (f + ke h + ne v ) C M 2 , e h = (3,0), e v = (0,V3). 
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Figure 1: Armchair graph for N = 10 and for N — 1. 

Vertices of T N are given by Tj + keh + ne v , (n, j, k) G Z x N 6 x Nat. If we identify the vertices 
?! + ne„ and ?x + iVe/j + ne v of for each nfZ, then we obtain the graph by 

= U^zT^, u = (n,j, k) G z = Z x N 6 x Zjv, Z w = Z/(JVZ), 

where = + ke h + we„, see Fig. [H [21 Let r w = r,- + A;e ft + ne v be a starting point of 
the edge T^. We have the coordinate x = + tej and the local coordinate t G [0, 1] on T u . 
Thus we give an orientation on the edge. For each function y on T N we define a function 
Vu = y\T U 7 u Z. We identify each function y w on with a function on [0, 1] by using the 
local coordinate t G [0, 1]. Define the Hilbert space L 2 (T N ) = ® ue zL 2 (T u ). Let C(T N ) be 
the space of continuous functions on r . We define the Sobolev space W 2 (T N ) that consists 
of all functions y = (y w ) we z G L 2 (T N ), (y^) u€ z G L 2 (T N ) and satisfy 
Kirchhoff Boundary Conditions: y G C(T N ) and for each vertex A ofT N satisfies: 

J2(-l)%(b)=0, where E A = {u G Z : A G T u }, b = b(u,A), (1.2) 

where if A = r w is a starting point ofT^ [i.e. t = a£ -A), i/ien 6(a;, A) = ; 
«/^4 = + ej is an endpoint ofT^ {i.e. t = 1 at A), then b(u,A) = 1. 
The Kirchhoff Conditions (11.21) mean that the sum of derivatives of y at each vertex of 
equals and the orientation of edges gives the sign ±. Our operator J4? on T N acts in the 
Hilbert space L 2 (T N ) and is given by (J^y)^ = -y" + qyw, where y = {y^^ez G £>(Jf ) = 
W 2 (T N ) and (y q y)uj = qyu, Q G L 2 (0, 1). Note that the orientation of edges is not important 
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Figure 2: A piece of a nanotube r . The fundamental domain is marked by a bold line. 

for the case of even potentials q e L 2 even (Q : l) = {q G L 2 (0, 1) : q(t) = q(l-t),t G (0,1)}. 
The standard arguments (see [KL] ) yield that J4? is self-adjoint. 

For the convenience of the reader we briefly describe the structure of carbon nanotubes, 
see |Haj . [S"DDj. Graphene is a single 2D layer of graphite forming a honeycomb lattice, see 
Fig. [31 A carbon nanotube is a honeycomb lattice "rolled up" into a cylinder, see Fig. [U In 
carbon nanotubes, the graphene sheet is "rolled up" in such a way that the so-called chiral 
vector Q = N\Q\ + N2Q2 becomes the circumference of the tube, where Q\, Q2 are defined in 
FigEl The chiral vector Q, which is usually denoted by the pair of integers (Ni, N 2 ), uniquely 
defines a particular tube. Tubes of type (N, 0) are called zigzag tubes. (N, A)-tubes are 
called armchair tubes. 

Our Schrodinger operator can be considered as a model Hamiltonain for 7r-electrons in 
armchair single-wall carbon nanotubes (SWCN) which attract tremendous interest from fun- 
damental science and technological perspectives [AL]. The best known model of SWCN is 
a discrete one obtained in the nearest-neighborhood tight-binding approximation proposed 
in 1992 (see the articles [HSU] . jMDWj . [SFDDj and review in [SDD] ) shortly after the dis- 
covery of carbon nanotubes [H]. Recently, the model was applied to SWCN (in particular 
to armchair configuration) subjected to an electric field [PBRLj . Being qualitative satis- 
factory, the tight binding approximation leads sometimes to quantitative errors, which are 
commonly addressed to curvature effects, see |OHKL| . However, as was shown in |RMTj . the 
tight-binding including up to third-nearest neighbors significantly improves approximation. 
To turn on interaction between all atoms in nanotube, one can proceed to the model con- 
sidered in |AT] , where carbon atoms was modelled by zero-range potentials on a cylinder. 
In the present article we use an intermediate model between the mentioned discrete and 
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Figure 3: The honeycomb lattice of nanotube. The unit cell is spanned by the vectors f^i 
and ^2- The type of the nanotube is defined by the pair (A^i,A^) € N 2 , N\ ^ N2, and 
corresponding chiral vector Q = N\£li + N 2 ^2- 

two-dimensional ones, namely we deal with quantum network models which attract lot of 
attention over the last decade, see |Kuj . The considered model was introduced by Pauling 
|Paj and was systematically developed in the series of articles by Ruedenberg and Scherr 
[RSJ. Further progress had been made toward periodic systems by Coulson in [Cj where a 
network model of graphite layer was worked out. A network model of a crystal with non- 
trivial potentials along bonds was studied by Montroll (Mj- Alexander [X] noticed simple 
relations between spectra of quantum graphs and combinatorial properties of discrete ones 
which became common for spectral analysis on quantum networks at the present time. 

The simplest pure carbon nanotube consists of atoms with covalent bond, when roughly 
speaking the potential of electric field between two atoms is even. The situation is more com- 
plicated, for example, in the model of a nitrogen-carbon-boron nanotube (see [Ha] . Ch.7.3), 
where the bond between atoms is ionic and the corresponding potential is non-even. 

Recall the needed properties of the Hill operator Hy = —y" + q(t)y on the real line 
with a periodic potential q(t + 1) = q(t),t G R. The spectrum of H is purely absolutely 
continuous and consists of intervals a n = [A^_ 1; X~],n ^ 1. These intervals are separated by 
the gaps 7 n = (A~, A+) of length \^ n \ ^ 0. If a gap 7 n is degenerate, i.e. \j n \ = 0, then the 
corresponding segments a n ,a n+ i merge. For the equation — y" + q{t)y = Xy on the real line 
we define the fundamental solutions A) and (p(t, A),tGl satisfying $(0, A) = <p r (0, A) = 
1,^(0, A) = <p(0,\) = 0. Let n = ^(1,0,0! = = = ^(1,0- The 

corresponding monodromy matrix A4, the Lyapunov function F, and the standard function 
F_ are given by 

The sequence Aq < A^ ^ A^ < ... is the spectrum of the equation — y" + qy = Xy with two 
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periodic boundary conditions, that is y(t + 2) = y(t),t G IR. Here equality A~ = A+ means 
that A^ is an eigenvalue of multiplicity 2. Note that F(X^) = (—1)™, n ^ 1. The lowest 
eigenvalue Xq is simple, -F(Aj~) = 1, and the corresponding eigenfunction has period 1. The 
eigenfunctions corresponding to A^ have period 1 if n is even, and they are anti-periodic, 
that is y(t + 1) = —y(t), t G M, if n is odd. The derivative of the Lyapunov function has a 
zero A„ in each interval [A~, A+], that is -F'(A n ) = 0. Let fi n , n ^ 1, be the spectrum of the 
problem — y" + qy = Xy, y(0) = y(l) = (the Dirichlet spectrum), and let v n , n ^ 0, be the 
spectrum of the problem — y" + qy = Xy, y'(0) = y'(l) = (the Neumann spectrum). Define 
the set <Jd = { u n,n ^ 1} and note that <Jd = {A G C : ip(l, A) = 0}. It is well-known that 
G [A n , A+], n > 1, and u ^ Aj. 
For simplicity we shall denote T a ^ C T 1 by T a , for a = (n, j) G Z\ = Z x Thus 
T 1 = U ae zi^a, see Fig [TJ In Theorem 11.11 we will show that M' is unitarily equivalent 
to H = ®i Hki where the operator Hk acts in the Hilbert space L 2 (r x ) and is given by 
{Hkf) a = - fa + qfa, {fa)aeZi, (fa)*eZi € ^(T 1 ), and the components f a ,a G Z x satisfy 
the Kirchhoff conditions: 

/n,l(l) = /n, 2 (0) = /n, 5 (0), /„, 2 (1) = /„, 3 (0) = / n>6 (0), 

Ml) = /M(0) = / n -i, 6 (l), s k j nA {\) = f nA (0) = / n _ li5 (l), s = e^, (1.4) 

/;,3W-/:,4(o)+ti, 6 (i) = o, s fc /:, 4 (i)-/;i(o)+/u5(i) = o. (1.5) 

Our analysis is close to previous papers [KLj, [KLlj about the Schrodinger operators 
on the zigzag graphs. We also introduce the monodromy matrix and the Lyapunov func- 
tions. We study the properties of these functions and here we essentially use the results and 
techniques from the papers |BBK] , |CK] , |K1] and |KL] . 

For the operator Hk we define fundamental solutions $ k = (v^ a ) a eZi, fi^ = (v 9 aTqW-Zh 
v = 1, 2 which satisfy 

— Ay + Yqy = Xy, on V , the Kirchhoff Boundary Conditions f ll.4p .( fl75|) 

Q k ,- X (l, A) = $'^(1, A) = I 2 , 6^(1, A) = $^(1, A) = 0, (1.6) 

where I n , n ^ 2 is the n x n identity matrix and 

Q k — I k,n,5 u k,n,5 j = I r k ,n,5 rk,n,5 j Q j\ 

\^k,n,6 ^k,n,6/ \Pk,n,6 fk,n,6/ 

We define the monodromy matrix (determined by the fundamental solutions on rj) by 

M k (X)=(^° !* )(1,A). (1.8) 
We formulate our first result about the monodromy matrix. 
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Figure 4: The supports of eigenfunctions: a) tp( 0,1 \ b) ^p ( -°' 2 \ 



Theorem 1.1. i) The operator is unitarily equivalent to H = ®±Hk- 

ii) Let k G Ztv- Then for any A G C \ od there exist unique fundamental solutions i?^ = 

(^Jaezi) <Pk = (PkJdaeZi, v=l,2. Moreover, each of the functions ^ a (x, A), (p^(x, A), 
x G T 1 is analytic in A G C \ op monodromy matrix .Mfc(A) satisfies 



detjM fc = l, Tr_M = 2(9F 2 - F 2 - 1), Tr jM fc = Tr A^ - 4ag, s fe = sin — , (1.9) 



and i/ie matrix-valued function TZM. k lZ 1 zs entire, where 71 = 7 2 © ^i-^- 

The monodromy matrix .A/U has poles at the points A G (Tb, which are eigenvalues of if^ 
(see Theorem II .21) . similar to the zigzag tube [KLj . [MVJ. 

Define the subspace 7ik(X) = {ip G Q(H k ) : iiZ&V' = A^>} for (A, A;) G cr pp (H k ) x Zjy. If 
dim7i fc (A ) = oo for some A G a pp (H k ), then we say that {A } is a flat band. In Theorem 
11.21 we describe all flat bands and corresponding eigenfunctions (see Fig. HI). 

Theorem 1.2. Let (X,k) G a D x Zjy. 27ien 

Every eigenfunction of7i k (\) vanishes at all verteces ofT 1 . 




(1.10) 




where 




(1.11) 
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ii) There exist two functions t/)^ '^ = (ipa )aeZ! £ ^fc(A),^ = 1,2, on T 1 snc/i i/ia£ 
supp C Uj =i r j ; each function ip^ 1 '^ = (i^n-m,j){m,j)eZi £ 7~^fc(A),n G Z. Moreover, 

each f G TCk(X) satisfies: 

f = ^(/n,l^ (n,1) + /n,2^ ( "' 2) ), (/n,l, /n )2 )n e Z 6 ^ © (1.12) 

ui/iere 

/n,i = Vi/;i(o), /n, 2 = /; 2 (o), «/ ^ 2 = i,fc = o, 

fn,l = : , fn,2 = : , *f<Px r 1 ° r <Pl = 1> k + 0, (1.13) 

H\H2 — V ><lX 2 — 1 

X] = 1 — s k ipf } x 2 = 1 — s^Vf. Moreover, the mapping f — > (/ n ,i, f n ,2)n& is a linear 
isomorphism between 7Yfc(A) and £ 2 © £ 2 . 

Using (11.111) we deduce that eigenvalues of A4k,k G Z^r have the form r^J, r^. 

Theorem 1.3. Lei t^t^ fre eigenvalues of Aik for some k G Zjy. Taen 
/%) Tne Lyapunov functions F k v = \(r k v ^ — — ), f = 1, 2 satisfy 

9^2 _ ^2 _ j 

Fv = ^-(-irV^> £fc = 4, Pk = (9F 2 -sl)4 + s 2 k Fl (1.14) 

Djfc(r, A) = det(M fc (A) - r/ 4 ) = (r 2 - 2F M (A)r + l)(r 2 - 2F K2 {\)t + 1), (1.15) 
where F^ v are branches of F k = ^ + ^/p^ on a two sheeted Riemann surface %Kf. defined by 
\fp~k- 

(ii) If Fk(X) G (—1, 1) for some A G K and A is not a branch point of Fk, then Fj.(X) ^ 0. 
(Hi) The following identities hold: 

a(H k ) = o-oo(H k ) U a ac (H k ), a^Hk) = a D , 

o-ac{H k ) = {Ag1: F kjV {\) G [-1, 1] for some v G N 2 }. (1.16) 

(iv) Ifq G L 2 even (0, 1), then a ac (H) = a ac (H ) = a{H). 

Definition 1. A zero of p k ,k G Zjy is called a resonance of H k . 

Let r kn) n ^ 1 be zeros of p k and let their labeling be given by Rer A T 1 ^ Rer^ $C 
R- er k~2 ^ ^ er k2 ^ ••• We snow the existence of real and non-real resonances. 

Proposition 1.4. Let q = q e = ^5(t — | — c k e — e 2 ),e ^ 0,t e [0, 1] for some k $ {0, y} ; 
and let no > 1. Then for each < n ^ there exist functions f kn {z 2 ), analytic in the disk 
{\z\ < e n } with some e n > 0, such that pk(r kn (£),q £ ) = for all e G (— e n , e n ) \ {0} and 

r± B (0) = (Trn) 2 , rg w (e) = (vm) 2 - 27m g ± ^^'V / 2 + 0(e 2 ) as v- 0. (1.17) 
^™er, w: n (e) = (rr n (e),r£ n (e)) C R and p fc (A,g £ ) < /or a// A G ct>r n (e), £ G (-e„,0). 




Figure 5: Graph of the Lyapunov function F k (X, q E ) for fixed e G (0, e-y). The spectral bands 
of H k (q £ ) are marked by bold lines, counted with multiplicity. Here £v>± are endpoints of the 
bands: E^ = r k2 and all other points -E^'± satisfy F k)U (E^) = (—1)", E^ n are periodic 
eigenvalues, E v ' 2n+1 are antiperiodic eigenvalues. 

Remark, (i) Graph of the function F k , k ^ {0, y} for q = q £ ,e G (0, e{) is given by Figj5l 

(ii) The potential q = q e = -5(t — \ — c k e — e 2 ), v ^ in Proposition 11.41 is not even. 

(iii) If e G (-e n , 0), then F k>v (X) <£ R for all A G uj^ n (e), which yields uj^ n (e) n cr(H k ) = 0. 
There are a lot of papers about the spectral analysis of the Schrodinger operator on 

periodic graph and periodic nanotubes. Molchanov and Vainberg |MVj consider Schrodinger 
operators with q = on so-called necklace periodic graphs. Korotyaev and Lobanov [KL], 
[KL1] consider the Schrodinger operator on the zigzag nanotube. The spectrum of this 
operator consists of an absolutely continuous part (intervals separated by gaps) plus an 
infinite number of eigenvalues with infinite multiplicity. They describe all eigenf unctions 
with the same eigenvalue. They define a Lyapunov function, which is analytic on some 
Riemann surface. On each sheet, the Lyapunov function has the same properties as in the 
scalar case, but it has branch points (resonances). They prove that all resonances are real 
and they determine the asymptotics of the periodic and anti-periodic spectrum and of the 
resonances at high energy. They show that there exist two types of gaps: i) stable gaps, 
where the endpoints are periodic and anti-periodic eigenvalues, ii) unstable (resonance) gaps, 
where the endpoints are resonances (i.e., real branch points of the Lyapunov function). They 
describe all finite gap potentials. They show that the mapping: potential — > all eigenvalues 
is a real analytic isomorphism for some class of potentials. 

Moreover, Korotyaev and Lobanov [KL1] consider magnetic Schrodinger operators on 
zigzag nanotubes. They describe how the spectrum depends on the magnetic field. 

Korotyaev [K2] considers integrated density of states and effective masses for zigzag 
nanotubes in magnetic fields. He obtains a priori estimates of gap lengths in terms of 
effective masses. 

Kuchment and Post [KuP] consider the case of the zigzag, armchair and achiral nanotubes 
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with q G Lg Uen (0, 1). They show that the spectrum of the Schrodinger operator (on these 
nanotubes), as a set, coincides with the spectrum of the Hill operator. We would like to note 
that the case of even potential q is simple and is closed to the case q = 0. Indeed, if q is even, 
then F_ = |MWj and by Theorem ll.31 the corresponding Lyapunov functions are expressed 
in terms of the Lyapunov function F (for the Hill operator). Moreover, using the identity 
F(X) = cosk(A) («(A) is the quasimomentum for the Hill operator), the Lyapunov function 
is expressed in terms of cosk(A), which is similar to the case q — 0, where F(X) = cos\/A. 

In contrast to [KuPj we provide spectral analysis for the Schrodinger operators with 
arbitrary potentials on edges, that can not be obtained using methods of |KuPj . 

We present the plan of the paper. In Sect. 2 we prove Theorem 1.2 about the eigenfunc- 
tions of Hk and Theorem 1.3 about a ac (Hk). A technical proof of Theorem 1.1 is placed to 
Sect. 3. 



2 Spectrum of the operators Hk 

Proof of Theorem 11.21 Proof of i) repeats the arguments from |KLj . 
ii) Recall x x = 1 — s k ipf } x 2 = 1 — s k ip± and ip[ ^ 0. 

If k ^ {0, y}, then Ims fc 7^ which gives Im j€\ ^ 0, Im x 2 ^ 0. 

If k = f , then s k = -1 and x x = 1 + <P? ^ 0, x 2 = 1 + <pf ^ 0. 

If k = 0, then s k = 1 and we have x x = 1 — ipf, x 2 = 1 — (pf. Thus, we obtain 

(a) If iff j£ 1 or k ^ 0, then x x ^ 0, x 2 ^ 0. 

(b) If <^f = 1 and k — 0, then x x = x 2 = 0. 

We will define the eigenfunctions of if^: in the case (a) ■0(°> 1 ) is given by 
4nf = 0, for all n^0,l,je N 6 , and = 0, j = 5, 6, ^ = 0, 4°f = x 2 <p, 

s k ^ff=^f= # = = -* Vi<C= -4°i 1} = = -*v? <C= «vsv, 

(2.1) 

and the function t// ' 2 ) is given by 

<f = 0, /orall (n,j)G(Z\{-l,0,l})xN 6 , ond ^ = 0, J = 1,3,4, 

<j 2) = 0, 4, <i 2) = 0, < 3 2) = **P, 

^% = -vas = s k ^f = # = -< 2 2) = ^& 2) = -^a a) = -^4T = m<p- 

(2.2) 

In the case (b) the functions ip(°' u \ v — 1, 2 are given by 

<f = 0, (n, j) G (Z \ {0, 1}) x N 6 , 4°f = 0, J G N 6 \ {4}, 

, (0,v) i (0,v) i (0,v) , i (0,1) , (0,1) / 
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4f = 0, < 4 1} = -^, <i 2) = < 3 2) = < 4 2) = 0, ^f = V . (2.3) 

Using (ED-©, we deduce that satisfy the Kirchhoff conditions (Ojl . fOjl . Thus 

?/;( ' r ') are eigenfunctions of if^. The operator iffc is periodic, then each ^ n,v \ (n, z/) G Z x N 2 
is an eigenfunction. We will show that the sequence ip( n ' u \ (n, z/) G Z x N 2 forms a basis for 
7ifc(A). The functions ?p( n ' u ^ are linearly independent: 

in the case (a) we have r„ j6 C supp?M n ' 1 ' ) \ supp and r nj3 C supp t//™' 2 ) \ supp ?/>( m ' 2 ) 

for all n 7^ m and r nj i C supp ip( n <^ \ supp ?/;( n ' 2 ) for all fiGZ; 

in the case (b) we have UfgNej^rnj C supp ip( n >"') \ supp ^ m ' u \ v — 1, 2, for all n ^ m, and 
r n 2 C supp^ ( ^ n ' 2 - ) \ supp^^' 1 -* for all n E Z,. 

Consider the case (a). For any / G Tlk(X) we will show the identity (11.121) . i.e., 

/ = /, where f = ^ /n> (n ^, and / n , 1; / n , 2 are given by (HH. (2.4) 

(n,i>)eZxN 2 

From A G a D , we deduce that f\ Vn] = f^(0)(p. The identity / = Y J fnA M and © 
provide 

/k.. = /n,l^ (n ' 1) |r„,e + /n,2^ (n ' 2) |r„, 6 = - /n,2)^, 

/k, 3 = /n.iV'^Irn,. + /n,2^ (n ' 2) |r„, 3 = (-/n,i + 
Substituting (11 . 131) into these identities we obtain 

fK 6 = fuM^ = /k, 6 , f\r n , 3 = f n Mv = /k, 3 a« n G Z. (2.5) 

This yields £ |/ n)1 ,| 2 < 00 and / G L 2 (r«). 

Note that / satisfies the Kirchhoff conditions ( 11. 4p . (11.51) and — + g/o, = A/ a , a G 
Consider the function u = f — f . The function u — at all vertices of T 1 and then 
u n,j = C n tj if, (n,j) G Z x N 6 , C n>3 = C n)6 = 0. The Kirchhoff boundary conditions (PUj)- 
ffOI) yield C n , 2 = C n , 4 = 0,n G Z. Assume that C nA = C. Then C n>b = ^ = g and 

Cn+1,1 = -^7^ = for all jigZ. Since it G L 2 (r 1 ), we have C = and it = 0, which yields 
(J23D- 

Consider the case (b). For any / G Hk(X) we will show the identity (11.121) . i.e., 

f = f, where f= £ / B , r ^ n ' r >, / n , 2 = /^(O), f nA = -^/^(O). (2.6) 

neZ,r=l,2 

From (12.61) and A G a D we deduce that 

/lr„,i = -/M^ = /|r nil , /|r n , 2 = /n,2<^ = /|r B , 2 all n el. (2.7) 

This yields l/n,^ 2 < 00 and then / G H(X). 

Consider the function u = f — f. The function u = at all vertices of T 1 , which yields 
u n,j = C n ,j¥>, n G Z, j G Z 6 , C n> i = C n;2 = 0. The Kirchhoff boundary conditions (jl.4B - fll.5p 
yield C n ^ = 0, and then C n ^ = 0, n G Z. Assume that C ni 3 = C. Then C n $ = — Co,3 = — C 
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and C n+ i^ = — C n .6 = C for all n E Z. Due to u E L 2 (T^), we have C = and u = 0, 
which yields ( 12. 6p . 

The mapping / — > (f n ,u)(n,u)eZxN 2 is a linear and one-to-one mapping from Hk onto 
£ 2 © £ 2 . Then it is a linear isomorphism. ■ 

Proof of Theorem 11.31 (i) Using the standard arguments from [KL] for Dk(-,r) = 
det(M k - rh) (see also [BBK] ) and (Ojl . (tUDj), OTj) we obtain (031) and (01) . 
The proof of (ii),(iii) repeats standard arguments from [KL]. 
(iv) If q is even, then F_ = (see p. 8, |MW] ) and identities ( 11.141) yield 

(3F + 1) 2 (3F-1) 2 , , 

*b,i = 1 2 ~ ~ 1, ^o, 2 = 1 2 J - 1, (2.8) 

Identities HE]) show that {Ael: F 0>V (X) E [-1, 1] for some v E N 2 } = {A E R : F(X) E 
[-1, 1]}, which together with fOf yields a ac (# ) = a(H). Identities (TETB1) . (TZ9D provide 

= {AGffi: F V (A) G [-1, 1] for some z/ E N 2 } C 

C{AGR: ^(A) E [-1, 1] for some z/ G N 2 } = ct„ c (# ) 
for all G Zat, which yields a ac (H ) = <j ac (H). * 

Proof of Proposition 11.41 If q E = ^5(t — a), e ^ 0, a E (0, 1), then we have 

sin zit — a) sinzt sinzasinz^ — a) 
v(t, A, q e ) = cos zt + cos za , (p(t, A, q e ) = h 



ez z ez 



s sin 2 „ sinz(2a — 1) i— 
F(X,q £ ) = cos z + —, F_(X,q £ ) = K — z = VX. 

Let a = \ + Cfc£ + e 2 , G" {0, y}. Then we obtain 

F_(A, &) = Sin2 " ( Cfc£ + g2) = c k + e + 0(e 2 ), F*_(X, q £ ) -c\ = 2c k e + 0(e 2 ) (2.10) 
2ze 

as |e| — > 0, uniformly on \z\ ^ 7TOq. Using the identities ( 11.141) . we rewrite the equation 
p fc (A, q e ) = 0, A G C in the form F 2 (A, g e ) - c 2 = and finally 

$ + (A,5)$_(A,e) = 0, g±(A, e) = g ( 3F(A ; fc 9e)Cfc ± V F - (A ' ge) ~ c ') ' AeC. (2.11) 

Using $±(A,0) = Ck 2 s s ^ z and ^±((7rn) 2 ,0) ^ 0,n > 1, and applying the Implicit Function 
Theorem to the functions $-t we conclude that for each n ^ n$ there exists e n > such that 
equation ( 12. lip has exactly two roots r^(e),r^(0) = (vrn) 2 , in the disc |e| < e n . 
Thus we have \fr±(e) = irn + t^, t = — > as |e| — > and then 

F(T^( g ),(fe) = (-1)" fcosr+ - 7 r 1 =^( £ + r + 0(r 2 )). (2.12) 



2(7rn + r)£/ 27me 
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Substituting (EZEDj) , (B3ZD into F 2 (\q £ ) -c\ = _ gg!^gsM ; we obtain 

2c fc £ + 0(e 2 ) = --^ 1 M , and then r = 0(e). (2.13) 

Using (jZJJ) again, we have r = -e ± gzEI f^ a v 73 ^ + 0(e 2 ), which yields f lTTTD . ■ 

3 Fundamental solutions 

Proof of Theorem 11.11 i) Define the operator S in by Su = (un,Ui, . . . ,wat_i) t , 
u = (u n )± G C^. The unitary operator S has the form S = s k V k , where Se k = s k e k 
and ek = — t-(1, s~ k , s~ 2k , s ~ kN+k ) is an eigenvector (recall s = e l ^~); VkU = e k (u,ek) is 

TV 2 

a projector. The function / in the Kirchhoff boundary conditions (11. 2\i is a vector function 
/ — (fu>),u = (n,j,k) G Z. We define a new vector- valued function f n j = (f n ,j,k)k=i ^ 
C N ,(n,j) 6 Zi = Z x N 6 , which satisfies the equation — + qf n j = A/ nj -, and the 
conditions for all n G Z, which follow from the Kirchhoff conditions (11.21) : 

/n,l(l) = /n, 2 (0) = / n , 5 (0), / n , 2 (l) = / n>3 (0) = / n>6 (0), 

/n,3(l) = /n,4(0)=/n-l,6(l), «S/ n , 4 (l) = /„,l(0) = / n _i, 5 (l), (3.1) 

4i(i) - /; 2 (o) - /; 5 (o) = o, /; 2 (i) - /; 3 (o) - /; 6 (o) = o, 

/i,3(l)-/;,4(0) + ti, 6 (l) = 0, 5^(1)-/; 1 (0) + /^ 1|5 (1) = 0. (3.2) 

The operators 5 and ^ commute, then we deduce that M"P k is unitarily equivalent to 
the operator H k acting in L 2 ^ 1 ) and H k is given by (H k f) a = —f£ + q(t)f a , where 
(fa)aeZ!, (fa)a£Z! G L^r 1 ) and components f a satisfy the boundary conditions (11 .4p . (11.51) . 
Thus M' is unitarily equivalent to the operator H = (BiHk 
ii) We prove the following idenities: 

(3.3) 

Let = y?(t, A), = <p'(t, A), ... Recall that any solution y of the equation — y" + qy = \y 
on [0, 1] satisfies 

»M = »<o)^<tfi)-*»<o)), M "=tt T)' (3 - 4 > 

Let .F n (i) = (/„, 5 (t), /n,eW, /n, 5 (0> /n,6(*)) T • We wil1 P rove tnat for each vector J"_i(l) = 
ft, = (h 5 , h s , h' 5 , h'o) 1 ^ G C 4 there exists a unique vector function (/ a )a6Zu satisfying the 
equations — + q(x)f a = \f a ,x G [0,1], a G Z ly and the Kirchhoff boundary conditions 

<HU), (O). 
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Firstly, we will determine /o,i,/o,3 in terms of h. The Kirchhof conditions fll.4p . (11. 5p 
at n = yield 

/o, 3 (l) = /o, 4 (0) = he, s k f 0A (l) = /o,i(0) = h 5 , 

/o, 3 (l) - + ^ = °> ^/W 1 ) " + fc B = 0- (3-5) 

Let Wl3 (t) = (/„,i(*), / , 3 (1 - 0, /o,i(*), /o,3(! " *)) T - Th en conditions Q imply 

u^b(0) = ( q 2 ?) ft + «; 4 , ^4 = (0, 0, ^, 4 (0)) T , jx = (J . (3.6) 

Using (JS3D, (JUS} we get / 0)4 (t) = Mt + ^{s~ k h 5 - $i/i 6 ), t £ [0, 1], which yields 

/o 4 (0) = — (s~ k h 5 - Ah), fU(l) = Mi + — (s-*ft 5 - ^ifte) = — (s-V^s - h 6 ), 

where we have used the identity "&i(p'i ~ $i<Pi = x - The last identities and (13. 6ft give 

^4 = (1 °) h, A = (j\ "^J , Wl3 (0) = F ft, Y = (| ?) . (3.7) 

Thus for each h = JF_ X (1) e C 4 there exists unique Wi 3 (0) = (/o,i(0), /o,3(l), /o,i(0), /o j3 ( i )) T 
and then the unique functions fo,i(t), /o,3(t),£ <G [0, 1] satisfying the equation —f" + qf = Xf 
and the Kirchhoff boundary conditions (13. 5ft . 

Identities « i m p ly (/$>) = * (jWQ , gg) = ^ g$>) , whieh 

yields 

^(d= W o), -j;), *■), j,=g 0.8) 

Secondly, we will determine / 0) 5, /o,6 m terms of / 0) i, /o, 3 - The Kirchhof conditions (11.41) . 
ffTol) at n = yield 

/o,i(l) = /o,a(0) = /o, 5 (0), /o, 2 (l) = /o, 3 (0) = /o >6 (0), 

fo,i(l) ~ /o l2 (0) - /o l5 (0) = 0, /i 9 (l) - /^(O) - /* 8 (0) = 0, (3.9) 

which gives 

^o(0) = ( q ?) wi 3 (l) + ^2, w 2 = (0, 0, -/i 2 (0), /^ 2 (1)) T , (3.10) 

where u; 13 (l) = (/o,i(l), /o, 3 (0), /^(l), /^(0)) T , ^o(0) = (/o, 5 (0), /o, 6 (0), /^(0), ^, 6 (0)) T . 
Using the Kirchhof conditions ([3]) and (13.41) . we obtain 

foAt) = ki(m + -(/o, 3 (0) - tfi/o,i(l)), £ 3 (0) = — (/o >3 (0) - 0i/o,i(l)), 
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/^(l) = / 0l i(l)^ + ^(/o,s(0) - ^i/o,i(l)) = -(^1/0,3(0) - / ,i(l)), 
where we have used the Wronskian $1^1 — fi'i^Pi = 1. Then 

W2=(i°a nW(l), A»=(*\ '})■ (3-11) 



±A 2 oj^» ^ 
Substituting this identity into (13.101) we obtain 

^o(O) = Y 2 u(l), Y 2 = (/^ . (3.12) 

Thus for each fteC 4 there exists a unique vector JF (0) = (/o,5(0), /o,6(0), /q 5 (0), /q 6 (0)) T 
and the unique functions /o,5(£), fofi(t),t <G [0, 1] satisfying the equation — /" + qf — Xf and 

the Kirchhoff boundary conditions ©. Identities give f{?' 5 JJ|) = -M f {°' 5 Jnll and 



v/;, 5 (o) 



fUW V/o, 8 (o) 



Thus we have proved that for each h £ C 4 there exist the unique functions foj(t),j G N 6 , 
satisfying the equation — /" + qf — Xf and the Kirchhoff boundary conditions (11.41) . ( 11.51) . 
By the periodicity, we obtain similar results for all functions f a , a E Z\. 

Recall that the fundamental solutions i?^ = (41)ae2u V 9 !^ = (rkDaeZn v — 1,2 satisfy 
(11.61) and the monodromy matrix M. k is given by (jl.8p . Then ^o(l) = M-kh, h = JF_ 1 (1), 
and identities lETTjl . (M . (|3~T2l) . <KW\ yield 

M fc = r 3 Win. (3.14) 

Now we calculate Y 2 YiY . We need the identities 

<p[ + i?i = 2F, tf'^i + Vi 2 = Vi^ - 1, + 0? = 2^iF - 1. (3.15) 

Identities (13T7L (JS3D give 

where we have used ( 13.151) and are given by (13.31) . Identities (I3.12p . (13.161) imply 
h 0\fV k <Pih\ = ( V k (pj 2 \ =v -i. 



(3.18) 
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where 



% ~ U»Vk + Ji^i A 2 +\b 2 ) ■ (3 - 19) 



Using <KTD\ . <KT6\\ . (EOTI) we have 

, T , . „ ( tix -l\ ( 2F -s k \ (\ 0\/2w'F-l -sVi i , , , 



-l 1 -ly v° 

Substituting the last two identities into (13.191) . we obtain the identity for T k in ( 13. 3ft . Substi- 
tuting identity (BTIBjl into (l3~T4j) . we obtain = Y z H- x T k ll which yields 1ET5|) . Identities 
(13.31) show that the matrix-valued function TZM.^RT 1 is entire. 

Identities ( 13.71) . (13.121) give det Y = det Y 2 = — 1. Identity (13.81) shows that Yi is obtained 

from the matrix by the transposition of rows and columns, hence detYi = 1. 

(I3.13P implies that Y 1 is obtained from the matrix Mr^j ^ ^ e ^ rans P os ^^ on °f rows 

and columns, hence det Y 3 = det.M det M.^ 1 = 1. Then (13. 14j) gives detTVU = 1. 
Furthermore, ( 13.31) gives 

Tr Vfc = Tr V = AF, Tr(V V k ) = 8F 2 + s~ k + s k . (3.20) 

Hence 



fej 



TrA^ fc = 4F($ 1 + ip' 1 )+8F 2 + S - k + s k -2 + 2ip 1 $' 1 = lQF 2 + s~ k + s k -2 + 2(p 1 $' 1 = Tr.M -4s 

where Tr M = 16F 2 + 2^i# / 1> which yields lOj) . 
We prove ( II. 101) . We have 

Tr.M 2 = Tr(tf 1 Vfc + V / VW 2 ) 2 + 2Tr(tf 1 / 2 + ^^^ 

(3.21) 

Direct calculations give 

Tr y fc 2 = Tr ^ 2 = 8F 2 + 2, Tr(V Vk) 2 = 32F 4 + 16F 2 (s k + s~ fc + 1) + s 2/c + s~ 2k , (3.22) 

Tr^xlfK) + <f'iV 2 V k ) = 32F 4 + 8F 2 (s fc + s~ k + 1). (3.23) 
Substituting identities (l3~20l . (|3~22l) and (l3~23l) into (l3~2Tj) we obtain 

Tr M 2 k =128F* + 8F 2 + 32F 2 tp 1 $' 1 + (32F 2 + 4^) {s k + s~ k ) + s 2k + S - 2k - 8^[ + 2^' 2 + 2 
= ^(Tr Mo) 2 + 72F 2 - 8s 2 k Tr M - As 2 2k - 4 = Tr M 2 - 8s 2 k Tr M - As 2 2k , 
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where we have used the identities s k + s k = 2 — 4sj*, #i</?i = ^i^'i + 1- This yields fll.lOp . 
Now we will prove (11.111) . Direct calculation show that Y T JY = J and 



where the identities j 2 Vfc = V fc T j 2 , h^o — Voh was used. Then (YT k ) T JYT k = J. Using the 
identities TZJTZ = tp 1 J, TZ^JTZ- 1 = ±J, we have ffTTTj) . ■ 
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